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ABSTRACT 

We report a sizable class of type 1 active galactic nuclei (AGNs) with unusually weak near-infrared 
(1-3/um) emission in the XMM-COSMOS type 1 AGN sample. The fraction of these "hot-dust-poor" 
AGNs increases with redshift from 6% at low redshift (z < 2) to 20% at moderate high redshift 
(2 < z < 3.5). There is no clear trend of the fraction with other parameters: bolometric luminosity, 
Eddington ratio, black hole mass and X-ray luminosity. The 3/xm emission relative to the 1/im 
emission is a factor of 2-4 smaller than the typical Elvis et al. AGN spectral energy distribution 
(SED), which indicates a 'torus' covering factor of 2%-29%, a factor of 3-40 smaller than required 
by unified models. The weak hot dust emission seems to expose an extension of the accretion disk 
continuum in some of the source SEDs. We estimate the outer edge of their accretion disks to lie at 
(0.3-2.0) x 10 4 Schwarzschild radii, ^10-23 times the gravitational stability radii. Formation scenarios 
for these sources are discussed. 

Subject headings: galaxies: evolution - quasars: general 



1. INTRODUCTION 

Characteristically, active galactic nuclei (AGNs) have 
hotter dust than starburst galaxies, which has long been 
employed to select AGNs in near-infrared (NIR) surveys 
(e.g., Miley et al. 1985, Lacy et al. 2004, 2007, Stern et 
al. 2005, Donley et al. 2008). In AGNs, dust reaches 
maximum temperature (^1000-1900 K) at the small- 
est equilibrium radius (Barvainis 1987, Suganuma et al. 
2006). The hot dust is assigned to the inner edge of the 
'torus' by the unified model of AGNs (Krolik & Begel- 
man 1988, Antonucci 1993, Urry & Padovani 1995). The 
obscuring torus might be a smooth continuation of the 
broad- line region (BLR; Elitzur & Ho 2009), with the 
BLR extending outward to the inner boundary of the 
dusty torus (Suganuma et al. 2006). 

However, exceptions have been observed and predicted. 
Observationally, Jiang et al. (2010) found two z ~ 6 
quasars, without any detectable emission from hot dust, 
suggesting that these hot-dust-free AGNs are the first 
generation quasars that do not have enough time to form 
a dusty torus. Rowan- Robinson et al. (2009) found 
AGNs with torus to bolometric luminosity ratio of only 
a few percent. 
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Emmering et al. (1992) proposed a BLR/torus struc- 
ture in the disk-wind scenario. The two structures 
correspond to different regions of a clumpy wind com- 
ing off the accretion disk rotating around the black 
hole (BH). The disk-wind scenario predicts that the 
torus disappears at luminosities lower than ~ 10 42 
erg s _1 because the accretion onto the central BH 
can no longer sustain the required cloud outflow rate 
(|Elitzur fe Shlosman 2006) . In a sample of nearby 
AGNs, even the BLR disappears at luminosities lower 
than 5 x 10 39 (M/10 7 Af o ) 2 / 3 erg s" 1 (jElitzur fc Ho I 
l2009lh 

The AGN structure is reflected in the shape of the 
spectral energy distribution (SED). The maximum dust 
temperature leads to a characteristic drop in the NIR 
emission at 1/im (Elvis et al. 1994, E94 hereafter; Glik- 
man et al. 2006). AGNs with weak or no 'torus' have an 
SED with weak or no IR bump, or no 1/zm inflection. 

We have found a substantial (~10%) population of 
such "hot-dust-poor" (HDP) AGNs by studying the 
SEDs of 408 XMM-COSMOS (Hasinger et al. 2007; Cap- 
pelluti et al. 2009) X-ray-selected type 1 (FWHM > 2000 
km s _1 ; Elvis et al. 2010) AGNs. All these X-ray sources 
have secure optical and infrared identifications (Brusa et 
al. 2010) and at least one optical spectrum, from either 
the Magellan (Trump et al. 2009a), Sloan Digital Sky 
Survey (SDSS; Schneider et al. 2007), or Very Large 
Telescope (VLT; Lilly et al. 2007, 2009) surveys. 

2. SELECTION OF HDP OBJECTS 

The HDP AGNs are selected based on the NIR and 
optical SED shapes. We plot NIR versus optical slopes 
on either side of lpm (rest frame, Figure 1; for details see 
Hao et al. 2010). The 1/itm point is not chosen arbitrar- 
ily. It is where the blackbody emission of the hot dust 
at the maximum sublimation temperature (e.g., 1500 K; 
Barvainis 1987) normally begins to outshine the emission 
of the accretion disk. 

The plot is equivalent to a color-color plot, but utilizes 
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more than four photometric bands. Briefly, we fit power 
laws (vF v oc v a ) on cither side of the 1/im inflection 
point of the rest-frame SED: l/im-3000A to derive an 
optical slope {chopt) and 3/im - 1/^m to derive a NIR 
slope {otNiR)- In normal type 1 AGNs, these slopes are 
a PT ~ 1 and a N iB ~ -0.7 (E94; Richards et al. 2006). 
In order to obtain slopes that reflect the shape of the SED 
correctly, we required the number of photometric points 
of the linear fit to be larger than 2. Four hundred four 
out of the 408 X-ray-selected type 1 AGNs 10 satisfy these 
criteria. The errors of the slopes are the standard errors 
of the linear fit. 

The aoPT—OLNiFt. plot can be conveniently used to sep- 
arate the nuclear from the host galaxy emission. Figure 1 
shows the aoPT and o.nir for the E94 radio-quiet AGN 
template (red cross) and three galaxy templates (a spiral 
— Spi4, a 5 Gyr old elliptical — E115, and a starburst — 
NGC 6090; big blue triangles) from the "SWIRE Tem- 
plate Library" (Polletta et al. 2007). The black curves 
show the slopes of the SED templates obtained by mix- 
ing the AGN and galaxy with different fraction (0%- 
100%) after normalizing the E94 and galaxy templates 
at 1/xm. The mixing curves of NGC 6090 and Spi4 define 
the boundaries of the possible slopes obtained by mixing 
the E94 with all the 16 available galaxy templates in the 
SWIRE library (small blue triangles). The magenta line 
shows the OPT and NIR slopes of E94 template when 
reddening {E{B — V) =0 — 1 mag) is applied 11 . Red- 
dening primarily affects aopr- 

The majority of the XMM-COSMOS AGNs (~ 90%) 
are explained by combinations of an AGN, a galaxy, and 
reddening (Hao et al. 2010), as they lie within the trian- 
gular region at the left of the otoPT—ctNiB, plot, bounded 
by the reddening and mixing curves. However, ~ 10% 
(41 sources) lie in the upper right corner of Figure 1, 
more than lcr beyond the AGN-host-reddening triangle. 
These sources have an optical slope consistent with typ- 
ical AGNs, but a relatively weak IR bump. We name 
them "hot-dust-poor" (HDP) quasars. 

3. HDP SEDS 

There is a range of SED shapes for the 41 HDP AGNs. 
We further divide these AGNs into three classes, accord- 
ing to their relative positions to the equal slope (red) line 
in the aopr — &nir plot. Figure 2 shows examples of 
SEDs for each class. All wavelengths in this section refer 
to the rest frame. 

Class I. Twenty-four sources lie below the equal slope 
line (blue symbols in Figure [lj . These sources have a 
normal big blue bump (BBB), but a factor of ~2-4 lower 
than the E94 mean SED at 3/xm. They all have a rather 
flat infrared SED shape (a nib. ~ 0) extending to at least 
10/im (rest frame; Figure [2j left). Severn have high lumi- 
nosities (> 5L*, ranging from 5.3L* to 24L„, where L* is 
the break point of the galaxy luminosity function; Cira- 
suolo et al. 2007) at 1/im and hence should have weak 
host galaxy contribution, while the rest have luminosities 

10 Three out of the other four (XID=2119, 5320, 5617) are low- 
redshift AGNs with only K band and one IRAC band in the rest- 
frame NIR SED. The remaining source (XID=54439) is the highest 
redshift AGN with only two IRAC bands in the rest-frame NIR 
SED. 

11 We applied the IDL dereddening routine ' 1 FMJJNRED . PRO ' 
with the SMC extinction curve (Gordon et al. 2003). 




_ g I i i i i I i i i i I i e_i i I i i i i I i i i i_ 

-3-2-10 1 2 

a OPT 

Fig. 1.— Slope-slope plot of the XMM-COSMOS type 1 AGNs. 
Red cross and red circle show the E94 mean SED and the lcr disper- 
sion of the E94 sources. The blue triangles show different SWIRE 
galaxy templates. The black lines connecting the E94 and the 
galaxy templates are mixing curves. The purple line shows the 
reddening vector of E94. The straight red solid line shows the 
olopt = &NIR line. Different colors of the points show different 
class of the HDP AGNs (I, II, and III, see the text for details). 
The black dots show all the other XMM-COSMOS type 1 AGNs. 



Class II. Six sources are consistent with lying on the 
equal slope line (green symbols in Figure[T]). The infrared 
emission could be the exposed continuation of the BBB 
to longer wavelengths (~2-3/^m). Two out of six, have 
high luminosities (> 5L*, 5.4L*-13L*) at 1/im, while the 
rest have luminosities 0. 322^-2. 7-L». 

Class III. Eleven sources lie above the equal slope line 
(red symbols in Figure [1]). These sources have flatter 
BBB than normal, possibly due to reddening. Some 
of them show quite strong ~ 10/im continuum emission 
(Figure [21 right). The class III sources have luminosities 
of 0.53L*-5.3L*. 

We calculated the mean SEDs of the HDP AGNs in 
each class and compared them with the mean SEDs of 
the other type 1 AGNs in the XMM-COSMOS sample 
having the same range of optical slopes for that class, 
i.e., normal type 1 AGNs with same BBB shape (Figure 
[21 upper). The HDP AGNs show relatively weaker 1- 
3/im emission, by a factor of 1.6 for class I, 3 for class II, 
and 2.5 for class III. 

4. SOURCE PROPERTIES 

4.1. Correlations 

We compared a number of observed properties of the 
HDP AGNs with the whole XMM-COSMOS sample: 
redshift (z), bolometric luminosity (L^oi, 24 /im- 912 A), 
BH mass (M BH ; Trump et al. 2009b; Merloni et al. 
2010), Eddington ratio (A#; Lusso et al. 2010), the op- 



Hot-Dust-Poor Quasars 




13.5 14 14.5 15 15.5 
logiz(Hz) (Rest) 



13.5 14 14.5 15 15.5 
logy(Hz) (Rest) 

MiixrfP. 1 01 



13.5 14 14.5 15 15.5 
logv(Hz) (Rest) 

1 0.1 




13.5 14 14.5 15 15.5 
logy(Hz) (Rest) 



13.5 14 14.5 15 15.5 
logy(Hz) (Rest) 



13.5 



14 14.5 15 
logw(Hz) (Rest) 



15.5 



TABLE 1 
K-S Probability 



Fig. 2. — Mean SED and examples of SEDs for each HDP classes: Upper row: mean SED of the three classes of HDP AGNs (black solid 
line) compared with E94 RQ mean SED (dotted red line) and the mean SED of the normal XMM-COSMOS type 1 AGNs with similar 
BBB (solid red line). Lower left: class I SED example (XID=2105, redshift z = 1.509); lower middle: class II SED example (XID=96, 
z = 2.117); lower right: class III SED example (XID=167, z = 2.048). The data points are from Capak et al. (2007), corrected for line 
emission as in Elvis et al. (2010). The SEDs are fitted with an accretion disk (red dashed line), a dust component (magenta dashed line, 
single temperature blackbody), and a galaxy component (blue solid line, 5 Gyr elliptical galaxy). The black dotted line shows the SED of 
the sum of the components. The black straight line shows the power-law fits that give appT an d &NIR- The host galaxy contribution is 
negligible, for XID=2105 and XID=96 at high luminosity (4L* and 13L* at 1/im). The green dashed line shows the blackbody fitting to 
outer edge of the accretion disk. 

tical to X-ray spectral index (a ox ) 12 , the X-ray lumi- 
nosity at rest frame 2keV (L^kev), the X-ray hardness 
ratio (HR) 13 , and radio loudness ((724) 14 • All the val- 
ues of these parameters were reported in Lusso et al. 
(2010), Brusa et al. (2010), and Hao et al. (2010). In 
the HDP sample, only one source 15 is radio-loud (Elvis 
et al. 2010). 

We used the Kolmogorov-Smirnov (K-S) test to com- 
pare the distributions of each parameter for the HDP 
AGNs with the whole sample (HDP AGNs excluded). 
The K-S probabilities are reported in Table 1. Only z 
shows a significant correlation (pk-s = 0.02%). Distri- 
butions of other parameters (Mbh, L2keV, logAs, and 
Lf, i) are indistinguishable for the two samples. We will 
study the X-ray properties of the HDP AGNs in a follow- 
ing paper. We can also see some of these results in the 
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where a ox < 0, and larger 



F„(2 keV) 
_F„(2500A)_ 

means X-ray louder (Tananbaum et al. 1979). 

13 HR= (H - S)/(H + S) where H are the XMM counts in the 
2-10 keV band and S those in the 0.5—2 keV band (Brusa et al. 
2007, 2010). 

14 <?24 =log(/24nm//i.4G.Hz) < as radio loud (Appleton et al. 
2004). 

15 XID=167, class III. 



fraction plot in Figure [3j where we divided two param- 
eters 16 into nine bins and plotted the fraction of HDP 
AGNs in each bin. 

The fraction of the HDP AGNs clearly increases with 
z. This fraction is 5.9% ± 0.2% at z < 2 but increases to 
20.4% ± 0.1% at 2 < z < 3.5. Beyond z = 3.5, the sam- 
ple is too small (one source in two bins) to give a useful 

16 One (z) with the lowest Pk—S anQl one (Lbol) with the highest 
PK-S' 
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Fig. 3. — Distribution in the redshift space (top) and fraction 
of the HDP in the whole sample as a function of redshift (middle) 
and Lf, ; (bottom). The black (number divided by 5) and blue 
histograms show the distribution of all the AGNs and the HDP, 
respectively. 

constraint. A linear regression fit of the fractions gives 
Ihdp = (0.031 ±0.002)(1 + )(i-2±o.i)_ The photometry 
used in the SED comes from a point source correction to 
the 3" aperture photometry. Low-redshift SEDs there- 
fore include less host compared to the high redshift ones. 
This effect will only move the source along the mixing 
curve, so it cannot explain the observed HDP fraction 
increase versus redshift. 

4.2. SED Fitting and Physical Parameters 

Studies of NIR spectrum of quasars showed that a 
blackbody spectrum is needed to fit the observed spec- 
trum (Glikman et al. 2006). The BBB is mainly caused 
by the emission of the accretion disk (e.g., Elvis et al. 
1994) and could be reddened by the dust in the host 
(e.g., Hopkins et al. 2004). At ~ l^m, host galaxy has 
the maximum contribution to the observed type 1 AGN 
SED (e.g., Polletta et al. 2007). Considering all these 
effects, we fit all the HDP AGNs observed NIR-optical 
SEDs with three components: (1) accretion disk emis- 
sion: we use the standard a-disk model in Schwarzschild 
geometry with electron scattering and Comptonization 
of soft photons in the disk atmosphere (Siemiginowska 
et al. 1995), with measured BH mass and accretion rate 
where available (Hao et al. 2010) and E(B-V) =0.1-0.2 
in eight cases to get good chi-square in the UV. (2)hot 
dust emission: we use a single temperature blackbody 
spectrum. (3)host galaxy emission: we use a 5 Gyr old 
elliptical galaxy (Polletta et al. 2007). Examples of the 
fit results are shown in the lower panel of Figure [2j 

4.2.1. Covering Factor 

The NIR SEDs are fitted with a blackbody to get the 
maximum dust temperature (Td). The sublimation tem- 
perature for graphite and silicate grains are ~ 1500-1900 
K and ~ 1000-1400 K, respectively (e.g., Salpeter 1977; 



Laor & Draine 1993). T d ranges from 800 K to 1900 K. 
Most (19 out of 34) HDP AGNs with dust component 
have Td >1400 K, suggesting graphite grains. 

Dust can only exist at a radius beyond the dust evapo- 
ration radius, given by r^ — l-3£„{ 2 46-^i50(? P c (Barvainis 

1987) , where L uv ,46 is the total ultraviolet (l/im-912A) 
luminosity in units of 10 46 erg s _1 and X1500 is the max- 
imum dust temperature in units of 1500 K. We find r^ 
ranges from 0.2 pc to 3.9 pc for the HDP sources. The 
emission area (Ad) is given from the blackbody fitting, 
and Ad ranges from 0.1 pc 2 to 10.68 pc 2 . The covering 
factor of the dust component is then f c = -Ad/(47rr?). 
For example, the results for source XLD=2105 (Figure[2J 
left) are 7^=1500 K and A<2=1.76 pc 2 . The evaporation 
radius is 0.83 pc, corresponding to f c = 20%. f c ranges 
from 1.9% to 29%, which are a factor of 2-40 smaller 
than the 75% expected from the unified models to give 
the observed type 2 to type 1 ratio (e.g., Krolik & Begel- 
man 1988). Even if considering the dependence of type 

2 to type 1 ratio on X-ray luminosity (e.g., Gilli et al. 
2007), the f c of HDP AGN is still smaller than the 50% 
expected for the X-ray bright AGNs. 

The above innermost radius estimation is compara- 
ble to those reported in Kishimoto et al. (2007, 2009). 
Kishimoto et al. (2007) also showed that the innermost 
radii measured by near-IR reverberation are systemati- 
cally smaller by a factor of ~ 3 than the value predicted 
with the Barvainis (1987) formula. For the HDP AGNs, 
if the innermost radius of the hot dust is really a factor of 

3 smaller, the covering factor would be a factor of 9 larger 
than the values we calculated, which could reach the 75% 
required by the unified models (e.g., Krolik & Begelman 

1988) for most (29 out of 34) of the HDP sources. 

4.2.2. Disk Outer Radius 

For 11 HDP AGNs with little host galaxy contamina- 
tion (defined as <50% at 1/^m), and available Mbh and 
\ E (Lusso et al. 2010; Trump et al. 2009b; Merloni et 
al. 2010), the lack of hot dust emission allows us to see 
what appears to be the accretion disk emission extending 
further out well beyond 1/mi. 

We fitted the outer edge of the accretion disk compo- 
nent (red dashed line in Figure [2]) with a single temper- 
ature (T c ) blackbody spectrum (green dashed line). The 
outer radius of the accretion disk (Frank et al. 2002): 

R out = 1.1 x 10 4 T~5 a-^r)-iM^\yi pc, 



1 { 6GAn 

1 ~ \TT?~i 



We 



where M 8 = M/(10 8 M Q ) and / = 

assume a = 0.1, r\ = 0.1. 

An example is the source XID=96 shown in Figure [2] 
(lower middle; M 8 = 3.55, Xe = 0.71), the blackbody 
temperature at outer edge is 3200 K. This gives R out — 
0.47 pc, ~ 1.3 x 10 4 Schwarzschild radius (r s ), 13.6 times 
the gravitational stability radius of the accretion disk 
(Goodman 2003, i.e., the radius beyond which the disk 
is unstable to self-gravity and should break up). For the 
11 HDP AGNs, the T c range from 2200 K to 4500 K. 
The R out range from 0.09 pc to 0.99 pc, corresponding 
to (0.29-2) x 10 4 r s . These R out are -10-23 times larger 
than the gravitational stability radius. 



5. DISCUSSION AND CONCLUSIONS 
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In a sample of 404 XMM-selected type 1 AGNs (ex- 
cluding four with incomplete NIR photometry) , we found 
that 10.1% ± 1.7% have weak NIR emission, indicating a 
relative paucity of hot dust emission. We call these HDP 
AGNs. 

HDP AGNs have not been reported from SDSS or ear- 
lier samples. We have made an initial check of where 
the Richards et al. (2006) and E94 samples lie on the 
otOPT ~ a NiR plane. We find similar fraction of HDP 
AGNs and will report on this study in detail in a later 
paper. The aopr — &NIR plane is a useful tool for lo- 
cating outliers. 

The fraction of HDP AGNs clearly increases to red- 
shift z = 3.5 (f HDP = (0.031 ± 0.002)(1 + z)( 1 - 2±0 - 1 )). 
No trends of the fraction of HDP AGNs with other pa- 
rameters, notably Li, i, are observed. 

We divided these HDP AGNs into three classes accord- 
ing to their atopr and o/vj.r- We fitted the HDP SEDs 
with three components: the accretion disk, a blackbody 
to represent the hot dust, and the host galaxy. We found 
the dust covering factors are 1.9%-29%, well below the 
typical 75% required by unified models (e.g., Krolik & 
Begelman 1988). For the 11 HDP AGNs with little host 
contribution (< 50% at 1 ^irn) and available Mbh and 
Xe, we estimated the lower limit of the accretion disk 
outer radius to be (0.29-2) x 10 4 r s (0.09-0.99 pc), corre- 
sponding to 10-23 times the gravitational stability radius 
of an a disk. How these disks stabilized is an open ques- 
tion, or the long wavelength turndown is not R out as we 
assumed. These results agree with the NIR disk spec- 
trum uncovered using polarized light in Kishimoto et al. 
(2008). We ignore further discussion of the sources with 
strong host galaxy contribution as the system is more 
complicated to include the also unknown galaxy part. 

There are several possibilities to explain the lack of 
NIR emission in HDP sources. First, most XMM 
COSMOS HDPs are at 1.5 < z < 3, when the universe 
is 2.1-4.2 Gyr old, and they have 1-3 Gyr from reioniza- 
tion to form a torus. The Jiang et al. (2010) proposed 
explanation thus seems unlikely for the XMM-COSMOS 
HDP AGNs. Second, the luminosity and Eddington ratio 
distributions of the HDP AGNs and the normal XMM- 
COSMOS type 1 AGNs are indistinguishable. This rules 
out the possibility that HDP AGNs are low-luminosity 
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or low-accretion rate sources, unable to support a dusty 
torus (Elitzur & Ho 2009). Third, evolutionary scenarios 
(e.g., Hopkins et al. 2008) predict short-lived Eddington 
limited outbursts after a merger, which could destroy the 
innermost dust. The HDP AGNs might then be quasars 
that had insufficient time to reform dust in the inner re- 
gion. Fourth, alternatively, when two galaxies merge , 
the chaotic dynamics might destroy the "torus". The 
'torus' then reforms in a timescale that might allow 10% 
of quasars to be HDP. The dust formation timescale de- 
pends strongly on pressure and temperature (Whittct 
2003; Krugel 2003) making estimates unreliable as both 
quantities are poorly known. Last, Guedes et al. (2010) 
suggested that a BH ejected by gravitational wave recoil 
and carrying along its accretion disk and broad emission 
line region, but not the hot dust "torus", would appear 
as an HDP AGN. 

The unusual SEDs of HDP AGNs can affect estimates 
of galaxy bulge and BH masses derived from SED fitting 
assuming the BH SED to be E94 mean SED. Merloni ct 
al (2010) recently used this method to show that z ~ 1.5 
quasars deviate from the local M—a relationship (Hiring 
& Rix 2004). Three of the Merloni et al. (2010) quasars 
belong to our HDP sample. For these sources, the SED 
fitting using E94 template gives a small AGN contribu- 
tion. Allowing an HDP SED to fit would increase the 
BH mass and decrease the stellar mass, making them 
deviate even more from the local M — a relationship. 
New AGN templates with variable dust bump strengths 
are needed to derive accurate galaxy and BH masses in 
these objects. 

The properties of the HDP AGNs still need to be in- 
vestigated to understand their formation scenario and 
put them in the context of galaxy and SMBH evolution. 
In particular, we will check the existing Hubble images 
and optical and X-ray spectra of HDP AGNs in future 
papers. 
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